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We investigate the properties of the one-dimensional Bose gas at zero temperature, for which exact 
results exist for some model systems. We treat the interactions between particles in the gas with an 
approximate form of the many-body T-matrix, and find a form of Gross-Pitaevskii equation valid in 
ID. The results presented agree with the exact models in both the weakly and strongly interacting 
limits, and interpolate smoothly between them. We also investigate the use of mean-field treatments 
of trapped BECs to describe the ID system in the strongly interacting limit. We find that the use 
of the many-body T-matrix to describe interactions leads to qualitative agreement with the exact 
models for some physical quantities. We indicate how the standard mean-field treatments need to 
be modified to extend and improve the agreement. 



PACS numbers: 03.75.Hh, 03.65.Nk, 05.30.Jp 



I. INTRODUCTION 

Recent advances in experimental techniques P, Q have 
stimulated interest in the topic of ultra-cold dilute Bose 
gases in one dimension. The one-dimensional (ID) case 
leads to some interesting physics which does not occur in 
higher dimensions, and is also of theoretical importance 
due to the existence of exactly soluble models 11,0,0- I n 
this paper we investigate the ID Bose gas in the limit of 
zero temperature. We extend our earlier work [|| on the 
interactions in two dimensions to the ID case, and obtain 
a ID form of Gross-Pitaevskii equation (GPE). We also 
apply standard methods of mean-held theory developed 
to describe three-dimensional Bose-Einstein condensates 
(BECs) to the ID trapped case. In both cases we com- 
pare the results obtained from these methods to the exact 
solutions. 

It is well known that a true BEC cannot occur at any fi- 
nite temperature in an interacting homogeneous Bose gas 
in ID 0, due to long wavelength fluctuations; nor can it 
occur in the limit of zero temperature |8j , due to quantum 
fluctuations. However, the presence of a trapping poten- 
tial changes the density of states at low energies, and in 
the weakly-interacting limit a BEC may be formed 0- 
In a recent paper, Petrov et al. discussed three different 
regimes of quantum degeneracy which can occur in a con- 
fined ID interacting system [lCj: BEC, quasi-condensate, 
and Tonks- Girardeau gas. In the weakly interacting limit 
a BEC exists, but as the interactions become stronger 
the mean- field energy becomes important. When this 
energy becomes of the order of the energy level spacing 
of the trap, fluctuations again become significant. In this 
regime, the system forms a quasi-condensate, with local 
phase coherence, rather than the global coherence asso- 
ciated with a true BEC. In this regime the gas is similar 
to the homogeneous case, in which it was shown that at 
T = the phase coherence has only a power-law decay, 



rather than the exponential decay occurring at higher 
temperatures [HI- Finally, in the limit of very strong 
interactions, the particles become impenetrable and the 
system is known as a Tonks- Girardeau gas 0,0. This 
limit is exactly soluble, due to a mapping relationship 
between the wave function for the bosons and that of a 
non-interacting fermion gas Q. 

In the following section we briefly outline the relevant 
exact results for a ID gas. Section lllll then presents 
a form of the GPE in which we use an approximate 
many-body T-matrix to describe the interactions in a 
ID gas. The GPE gives the correct results in both the 
weakly-interacting and Tonks-Girardeau limits, and in- 
terpolates smoothly between these regimes. We confirm 
that the non-linear dependence in the GPE changes from 
cubic in the weakly-interacting regime to quintic in the 
strongly-interacting regime, and that this arises due to 
the many-body effects on the scattering of particles. In 
Sec. II VI we then discuss several standard mean-field the- 
ories and present the results of these for the trapped ID 
gas. In the weakly interacting limit these theories should 
be valid since a true BEC exists in this regime. Our aim 
is to push several different theories into the strongly- 
interacting limit with the dual objectives of describing 
the ID system and investigating the strengths and weak- 
nesses of these theories. 



II. EXACT RESULTS IN ID 

The ID homogeneous system of N bosons which inter- 
act via a delta-function potential of the form 



V(x\ - sc 2 ) = V 5(xi - x 2 ), 



(1) 



is exactly integrable. In the limit Vo — > oo the particles 
are impenetrable (the Tonks-Girardeau gas), and exact 
results can be obtained even in spatially inhomogeneous 
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systems. The exact solutions are found by applying the 
appropriate boundary conditions to the wave function 
describing the system. 

In the impenetrable limit the wave function must nec- 
essarily vanish wherever two particles meet (i.e. wher- 
ever Xi = Xj). Girardeau showed in [3j that the exact 
A^-boson ground state wave function subject to such a 
boundary condition is equal to \ipF(xi,...,xjf)\, where 
i/jf(xi, ...,Xn) is the wave function of iV non-interacting 
spinlcss fermions governed by the single-particle Hamil- 
tonian H sp (including terms from any external potential). 
Since the wave functions of such fermions must be anti- 
symmetric with respect to co-ordinate exchanges, their 
wave functions automatically vanish whenever xi = Xj, 
and so the Tonks-Girardeau gas boundary condition is 
implicitly obeyed by such a gas. 

A number of bulk properties of a ID Tonks-Girardeau 
gas may be calculated quite trivially as a result of this 
Fermi-Bose mapping theorem |3j, Il3| . For example, the 
ground state energy of a Tonks-Girardeau gas is sim- 
ply found from a summation of the energies of the first 
N non-interacting single-particle states (i.e. those that 
would be occupied by N ideal fermions in the ground 
state configuration). Thus, for a homogeneous gas Q 
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N 
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m b 



(2) 



where n is the density. The chemical potential is then 
easily calculated as 



dN 



m 2 ' 



(3) 



Similarly, the chemical potential for a Tonks-Girardeau 
gas in a harmonic potential of frequency u) is 



fi = huj(N + 1/2). 



(4) 



The system of N bosons interacting with an arbitrary 
value of V was solved in 1963 by Lieb and Liniger 0, || 
for a homogeneous gas. The Tonks-Girardeau gas results 
are recovered in the limit Vo — > oo. The ground-state 
energy of the Lieb-Liniger gas is given (at T = 0) by 



Eq ft 2 

~N = 2m" 6(7) ' 



(5) 



where e("f) is a solution to the Lieb-Liniger system of 
equations, and is a function of the dimensionless quan- 
tity 7 = Vom/h 2 n which was found to be the important 
parameter in the system. The function e("f) has the fol- 
lowing asymptotic behaviour 



e{i) 



7 

nl ( i \ 

3 



for 7 
for 7 



0. 

oo, 



(6) 



and the energy per particle agrees with Eq. in the 
impenetrable limit. The chemical potential is given by 
the equation 



2m 



3e(7) - 7 



dej-y) 
c?7 



(7) 



In the 7^0 weakly-interacting limit, the exact results 
were found to be in agreement with the predictions of 
standard Bogoliubov theory. The Bogoliubov predictions 
begin to fail as 7 is increased and for t> 2 the results 
do not agree well with the exact theory |j] . Interestingly, 
for fixed Vq the 7 — > limit corresponds to the high- 
density limit, whilst a high value of 7 implies a small 
density. This is then the reverse of the 3D case where 
Bogoliubov theory is valid in the low density limit. In a 
subsequent section we will show to what extent higher- 
order modifications to Bogoliubov theory can improve 
the agreement for higher 7. 

For a system where we keep the density constant, 
and increase the interaction strength Vb from the non- 
interacting limit to the impenetrable limit (i.e. increas- 
ing 7), we expect the chemical potential to vary in the 
following manner. In the low-Vo limit the chemical po- 
tential increases linearly with Vo, but as Vo is increased 
further fi starts to saturate, and in the high-Vo limit it 
becomes independent of Vb, in agreement with Eq. J3J. 
In the following section we introduce a modified GPE 
which follows this behaviour. 



III. THE GPE IN ID 

We consider the model system of a gas of one- 
dimensional point particles interacting with an arbitrary 
strength Vo, corresponding to the Lieb-Liniger gas. We 
emphasize here that we are interested in the limit in 
which the interparticle interactions in the system are gen- 
uinely one-dimensional, because it is only in this limit 
that we can compare to the exact results mentioned ear- 
lier. 

We can write a one-dimensional Gross-Pitaevskii equa- 
tion for a system of N particles by a straightforward gen- 
eralisation of the usual 3D equation, giving 



n2 d2 ,r , x 
"2r^ + Wa;) 



-g m N\^{xt 



ip{x), 



(8) 

where ip{x) is the condensate wave function, and Vt rap (a;) 
represents the trapping potential. The coupling parame- 
ter (/id describes the interactions between particles, and 
we discuss the nature of this quantity in the following sec- 
tion before presenting the results predicted by the GPE. 



A. Interactions and the GPE 

The GPE can easily be derived by functional differ- 
entiation of the Hamiltonian for a gas with a general 
interparticle interaction V(x). The coupling parameter 
is then given by the matrix element (0|V(x)|0), where 
the particles in both the incoming and outgoing states 
have zero relative momentum (see, for example, [14|). 
However, a representation in terms of matrix elements 
of the bare interparticle potential is undesirable, since 
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for a realistic potential they are often very large and fre- 
quently the exact potential itself is unknown. In order 
to avoid this it is common for 3D problems to use the 
two-body T-matrix T 2 b in place of the interparticle po- 
tential. This T-matrix describes the scattering of two 
particles in a vacuum, and it is possible to rewrite the 
interaction term in the Hamiltonian in terms of T 2 b |14| . 
Carrying out the functional differentiation again leads to 
a GPE of the above form, where the coupling parameter 
is now g = (0\T 2B {E = 0)|0). Here E is the energy of 
the interacting particles, and for two condensate particles 
this is zero. In three dimensions this leads to the familiar 
result (?3D = A?:h 2 a/m, where a is the s-wave scattering 
length. However, in ID (and 2D) the two-body T-matrix 
vanishes in the zero energy and momentum limit, and has 
significant imaginary terms at a finite real energy (see, 
for example, [la]). 

In earlier papers [fillip, we have argued that the many- 
body effects of the surrounding atoms on particle interac- 
tions, while being a small perturbation in 3D, are crucial 
in lower dimensions. The many-body effects can be incor- 
porated into a many-body T-matrix Tmb [3 EH 
and the coupling parameter in the GPE becomes 



9iD = (0|T MB (£ = 0)|0). 



(9) 



Unfortunately, the many-body T-matrix is very difficult 
to solve exactly. However, in |6J] we presented an ap- 
proximate solution, valid in the zero temperature limit, 
in which the many-body T-matrix was approximated by 
the two-body T-matrix evaluated at a shifted effective 
interaction energy (off the energy shell). The GPE in ID 
is therefore given by Eq. JSJ with the coupling parameter 

9m = (0\T M b{E = 0)|0) « (0|T 2B (£*)|0), (10) 

where E* is the shifted interaction energy, and was found 
in [6| to be negative. As an illustration of the utility and 
accuracy of this approximation, we discuss its application 
to the 3D case in the appendix. 



B. The hard-core limit 

We will look first at the GPE in the impenetrable, 
Vq — > 00, limit in which the system is a Tonks-Girardeau 
gas. The off-shell two-body T-matrix for such a ID gas, 
calculated in 0, is 

(0\T^ ks (E)\Q) = 




(11) 



where a is the size of the impenetrable particle, and for 
the delta- function gas we take the limit that a — > 0. Com- 
bined with Eq. (|10l) this gives a coupling parameter of the 
form M 



ffTbnks(A<) = \/ 4C?i ji/rn 



to leading order for a homogeneous system. Here we 
have taken E* — —Cfi and assumed that C > 6], (note 
that if this is not the case then the coupling parameter 
is imaginary to leading order). Using the homogeneous 
solution to the GPE /i = ^Tonks^ we can rewrite the 
coupling parameter in terms of the density n as 



ACh 2 
STonks(n) = n. 



(13) 



Combining this with fx = .gTonks^ gives \i oc n 2 , and a 
comparison with the exact result found by Girardeau and 
given in Eq. J2J), shows that we have obtained the correct 
functional dependence of fi on n. This stems from the 
dependence of the coupling parameter on yfjj, in Eq. I|12(l . 
This is further verification of our approximation for the 
many-body T-matrix in terms of the off-shell two-body 
T-matrix. A comparison with the Girardeau result also 
shows that C = 7r 2 /8. 

To apply this result to a trapped system, we use a 
local density approximation, assuming that the scatter- 
ing between two particles occurs on a length scale much 
shorter than the variation in the potential. First we solve 
the homogeneous expression [i = (?Tonks(A t ) rl to obtain an 
expression for the chemical potential in terms of the den- 
sity n. Applying the local density approximation then 
requires /x(n(x)) = gTonks(x)n(x) to hold for all x, and 
we obtain a spatially-dependent coupling parameter 



. TT 2 h 2 N. 2 
ffTonksW = — \V{X)\ , 

2 m 



(14) 



where we have used n(x) = N\ip(x)\ 2 . 

We can now write the GPE for a ID trapped conden- 
sate as 



jj,ip(x) 



h 2 d 2 . 

^^ + W*) + 



TT 2 h 2 N 2 

2m 



ij)(x). 
(15) 

This GPE differs significantly from the 3D version by 
virtue of the fact that the nonlinear term is now propor- 
tional to \ip\ 4 ip instead of |?/'| 2 ?/'- In the Thomas-Fermi 
limit we can neglect the kinetic energy term and the so- 
lution for the wave function can easily be shown to be 



\^{x)\ 



1/4 



(P - Utrap(z)) 

{h 2 ir 2 N 2 /2my/ 4 



#(/•«- U trap (a;)) 



(16) 



(12) 



This differs substantially from the 3D and 2D cases, 
where the Thomas-Fermi wave function is instead pro- 
portional to (/j — V^ap^)) 1 / 2 , due to the different form 
of the non-linearity. 

The Tonks-Girardeau GPE and the results presented 
in this section are in exac t ag reement with the earlier 
work of Kolomeisky et al. |20|, who obtained the same 
form of coupling parameter from a renormalisation group 
approach [21||37l|. 

In Fig.QJi we present a comparison of the density iV|V>| 2 
obtained from the solution of the GPE of Eq. I|15|l with 
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many-body T-matrix in a homogeneous ID system is 



5 

x (units of ^trap) 



FIG. 1: (a) Comparison of ground state density for Tonks 
gas (solid) and density calculated from Eq. 115H with |'i/ , | 4 V' 
non-linearity (dashed), both with /i = 20hu and N = 20. 



the exact result from the Bose-Fermi mapping theorem 
calculated by a summation of the first N single-particle 
states of the trap. It can be seen that the GPE solu- 
tion is in very good agreement with the exact result, and 
differs only in the absence of the oscillations seen in the 
exact result. These oscillations rapidly average out as 
the number of atoms is increased, and for a moderately 
large number the difference between the GPE results and 
the exact results will be negligible. 



C. The delta-function potential 

We now turn to a discussion of the Lieb-Liniger gas. 
We are interested here in describing the gas from the non- 
interacting limit right through to the Tonks-Girardeau 
limit by means of an appropriate GPE. 

The same argument which we used in deriving the 
Tonks-Girardeau GPE can again be used, namely that 
the coupling parameter should be the many-body T- 
matrix. Again we approximate this by a two-body T- 
matrix evaluated off the energy shell at an energy of 
— Cfi, where we use the value of C — ir 2 /8 obtained 
above. 

The two-body T-matrix for a delta-function potential 
is given by 0] 



<0|T 2 ^(£)|0> = 



V 




(17) 



Note that for Vq — * oo, this reduces to the Tonks limit 
of Eq. (|11|) . It follows that our approximation for the 







1 + Vb 



2m 1 



(18) 



This energy-dependent form of the coupling parameter 
can be transformed to a density-dependent expression 
by means of the local density approximation, as outlined 
above. This leads to the expression 



1 



(19) 



n{x) 



V 2 + V n(x) - 



2ir 2 n 2 



V 



which describes an inhomogeneous system. Note that 
this expression gives the correct results for the homoge- 
neous gas in the extreme limits Vq — * and Vq — » oo, 
giving /i = nVo and Eq. respectively. 

Thus we can now use the above density-dependent ex- 
pression for the coupling parameter in Eq. © to obtain 
a GPE for a trapped ID gas for any arbitrary Vo- 

We use the GPE here to describe the system over the 
full range of Vq, despite the fact that, as mentioned ear- 
lier, a true BEC only exists in a trap in the weakly- 
interacting limit pH |. Once the interactions cause suf- 
ficient perturbation to the density of states that their 
discrete nature is smeared out then no condensate will 
form in ID, and in the very high-Vo (Tonks) limit it is 
known that there is no condensate. We can still use the 
GPE in the high-Vb regime however, as our previous re- 
sults have shown, provided that we now reinterpret the 
"wave function" as simply the square root of the density 
ifj(x) — y/n{x), rather than the coherent part of the field 
operator. 



1. The Homogeneous Case 

In the homogeneous case we can compare the results of 
the GPE with the interaction strength given by Eq. I|19|) 
to the exact results for a Lieb-Liniger gas. In Fig. we 
compare the chemical potentials resulting from the GPE 
with the exact result of Eq. J3J). The results from the 
GPE show the correct behaviour in the two extreme lim- 
its of 7, and the overall agreement is quite good. How- 
ever, the chemical potential of the exact results can be 
seen to rise more steeply in the intermediate region. We 
note that the results have been obtained using a value 
of C which strictly applies only in the large 7 limit. A 
larger value of C gives a curve which agrees very well 
with the exact results in the weakly-interacting regime, 
but predicts the wrong value in the large 7 limit. It seems 
possible therefore that the value of C should actually be 
dependent on 7, but a detailed investigation of this is 
outside the scope of this paper. 
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FIG. 2: The chemical potential as a function of the parameter 
7 for the homogeneous Lieb-Liniger gas. The dashed line 
shows the exact solution from the Lieb-Liniger equations Q, 
while the solid line shows the prediction of the GPE with 
the coupling parameter given in Eq. 118H . Note that in the 
extremely high 7 limit (not shown) both lines have the same 
asymptotic value. 



FIG. 3: Solutions to the GPE for a ID gas of 15 atoms inter- 
acting via a delta-function potential. The dashed curve corre- 
sponds to the weakly-interacting limit with Vb = 0.2ftw£trap. 
The dot-dash and dotted curves correspond to Vb = 3?iu>£trap 
and Vo = 17hu}£tia,p respectively, while the solid line is close to 
the impenetrable gas limit with Vb = 250hui£ t Ta, P - These val- 
ues correspond to 7t rap = 0.02, 0.63, 5.3 and 99 respectively. 



2. The Trapped Case 

We now turn to the case of a trapped gas, where 
we will parameterise our results both by the interac- 
tion strength Vo and by the dimensionless parameter 
7trap = VQm/h 2 no where no is the peak density of the 
system (analogous to the parameter used by Licb and 
Liniger). 

We expect that in the low-Vb (small 7trap) limit the 
chemical potential will vary as 

u«V 2/3 , (20) 

which is the Thomas-Fermi result obtained from the GPE 
with a cubic non-linearity (this is not quite the cor- 
rect prediction in the extremely low-Vb limit where the 
Thomas-Fermi approximation is invalid, but it should de- 
scribe the majority of the low-Vb regime adequately) . In 
the opposite high-Vb limit, the Bose-Fermi mapping the- 
orem for the trapped Tonks-Girardeau gas led to Eq. Pfl 
which is completely independent of Vb- Note that our 
derivation of the GPE in Eq. (JSJ neglected terms of or- 
der 1/iVo since it was derived in the high-iVo limit. For 
this reason we should expect only to reach an asymptotic 
chemical potential of fi = NqHuj, which is approximately 
equal to Eq. (0} in this limit. 

Figure[3]shows sample solutions to the GPE for a range 
of Vb for iVo = 15. In the weakly-interacting limit the 
wave function can clearly be seen to be a gaussian, as 
expected. However as the interactions are increased the 
wave function broadens until in the high-Vb limit we have 
the shape discussed earlier in the hard-core case. Fig- 
urelUshows the increase in the chemical potential over the 
same range. As expected from the exact results for the 
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50 100 150 200 250 



Vb (units Of TlCi^trap) 

FIG. 4: The chemical potential for a gas of 15 atoms interact- 
ing via a delta-function potential as a function of the potential 
strength (solid line) and the parameter 7t rap (dashed line and 
upper axis). 

hard-core gas, in the high-Vb limit the chemical potential 
becomes almost independent of the interaction strength. 
For the data shown the chemical potential does not ap- 
pear to quite reach the limit pt = NqHuj expected from 
the exact results. However, in a calculation with the 
extreme value of Vb = 10 4 fi.w^trap we obtained the result 
(i = 15.00?iu; as expected. This agrees with the behaviour 
found by Lieb and Liniger in the homogeneous case, who 
found that the results reached the asymptotic limit very 
slowly (e.g. for 7 ss 4 the chemical potential had only 
reached ~ 0.9 of its asymptotic value 3). 

We have therefore presented a form of the GPE in 
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which the coupling parameter is an approximate many- 
body T-matrix given by the off-shell two-body T-matrix. 
We have shown that it agrees with the exact results 
for the Lieb-Liniger gas in the weakly-interacting and 
strongly- interacting limits, and it provides an interpo- 
lation between these limits in the intermediate regime. 
A drawback of using this GPE is that we lose any in- 
formation about the coherences involved in the system, 
since we reinterpret the equation simply in terms of the 
density in the high-Vo limit. The concept of using the 
many-body T-matrix in the GPE stemmed initially from 
higher order theories in which it appears naturally via the 
anomalous average (defined below). We therefore turn, 
in section IIVI to an exploration of these theories in ID 
and see how the changes in the interaction strength arise 
naturally. 



3. Quasi-ID gases 

Although we have been dealing with the case of a 
model delta-function potential gas, the results can to 
some extent be mapped onto those for a quasi-lD system. 
In practice, a ID system could be created by squeezing 
a 3D trapped gas condensate in two dimensions. As the 
confinement is increased, the system progresses from be- 
ing 3D to being ID, via an intermediate regime we call 
quasi-lD for the purposes of this paper. In this interme- 
diate regime, the dynamics of the system as a whole are 
frozen out in two dimensions, but the scattering between 
particles is not yet genuinely one-dimensional. 

In the 3D regime the system is described by a GPE 
with a IV'I 2 ^ nonlinearity, as for the low Vq Lieb-Liniger 
gas. However, in the strict ID limit the atoms will be 
unable to scatter past one another, and will form a Tonks- 
Girardeau gas. In this limit the appropriate GPE will be 
that of Eq. (|15fl . The progression from a 3D gas to a ID 
gas therefore mirrors in many respects the progression 
which occurs if the potential strength Vb is increased for 
the Lieb-Liniger gas. The 3D to ID progression has been 
investigated in more detail in references [2^, |2^, 0] . 



IV. HIGHER ORDER THEORIES IN ID 

In the previous section we used the GPE in conjunc- 
tion with an approximate many-body T-matrix in order 
to describe the Lieb-Liniger gas. In this section we ex- 
amine the suitability of higher-order theories to describe 
this gas. Ideally we would like a theory which makes the 
transition between the weakly interacting and strongly- 
interacting regimes smoothly, and which accurately de- 
scribes the coherence of the system as well as the den- 
sities. There are several higher order theories which are 
used to good effect to describe 3D condensates, and in 
this section we will apply some of these to the ID trapped 
gas case. 



Since we can always choose a value of Vq to be suffi- 
ciently weak that a condensate does exist in a trapped 
zero-temperature ID gas, there will always be a regime 
in which the 3D theories can be applied. As we increase 
Vq the condensate is replaced by a quasicondensate with 
only short range phase coherence, and this will lead to 
problems in the higher order theories. However, we have 
described in the previous section the behaviour that we 
expect in the high-Vb limit. It is therefore of interest to 
apply the higher order theories to see how and why they 
fail in this limit, and to gain insight into any improve- 
ments needed to describe the ID gas adequately. 

As we have seen, the many-body effects on interac- 
tions are fundamental in understanding the behaviour of 
low dimensional BECs. We must therefore use a the- 
ory in which these effects are accounted for accurately. 
Such a theory is provided b y th e gapless-Hartree-Fock- 
Bogoliubov (GHFB) theory 1!>. 25]. In this theory, the 
condensate wave function if>o(x) is described by the gen- 
eralised GPE in ID 

h 2 d 2 

Hipo(x) = ~y~ ^2^o(r) + V r tra P (a;)'0o(a;) 
+ N U con (x)\<ipo(x)\ 2 ip (x) 
+ 2U cx (x)p(x)Mx)- (21) 

Here p(x) is the non-condensate density and it is defined, 
along with the associated anomalous average k [needed 
in Eq. by 

p{x) = £(K(z)| 2 + | Vi (:r)|Vi + Mz)| 2 , (22) 

i=£0 

K (x) = 5^u i (xK(a;)(2n i + l). (23) 

i=£0 

The Bogoliubov functions Ui(x) and Vi(x) are the solu- 
tions of 

C(x)ui(x) + M{x)vi(x) = EiuAx), 

(24) 

C(x)Vi(x) + M*(x)ui(x) = -£<«i(x), 
solved in a basis orthogonal to ipo(x), and with 

L(x) = H sp - p + 2N U con (x)\Mx)\ 2 

+ 2U ex (x)p(x), (25) 
M(x) = N Q U con (x)i; 2 (x). (26) 

The quasiparticle population factors rii are given by the 
Bose-Einstein distribution function, and for the zero tem- 
perature case we are concerned with here these factors 
vanish. 

This theory distinguishes between the interactions 
which occur between two condensate atoms (governed 
by U con (x)), and those which occur between a conden- 
sate atom and a non-condensate atom (governed by U ex ) . 
The various higher order theories which we will explore 
in this section can now be summarised thus: 
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BdG By setting U cx = (i.e. neglecting interactions be- 
tween the condensate and the non-condensate), and 
using U con (x) = Vq, we obtain the familiar GPE 
and Bogoliubov-de Gennes equations. 

HFB-Popov The HFB-Popov theory corresponds to 
the choice U con (r) = Vq and U ex (r) = Vo. 
Thus both condensate-condensate interactions and 
condensate-non-condensate interactions are in- 
cluded, although many-body effects on the inter- 
action strength are not. 

GHFB1 The first version of the GHFB theory is given 
by setting U con (r) = T MB (r) [defined in Eq. (|2?I) ] 
and U cx (v) = Vq. This version of GHFB therefore 
improves upon the description of the BdG or HFB- 
Popov descriptions by including many-body effects 
on the scattering of two condensate atoms. 

GHFB2 From the above arguments it is clear that 
another GHFB theory may be proposed. This 
consists of setting U con (r) = Tmb(i*) and 
U ex ( r ) = r MB (r), which upgrades the descrip- 
tion of both condensate-condensate interactions 
and condensate-non-condensate interactions to the 
many-body T-matrix. 

In the GHFB theories, the many-body T-matrix is 
given by [lid 

Tmb(x) =Vo(l + -^-), (27) 

and so we can see that the inclusion of the anomalous 
average leads to the inclusion of the many-body T-matrix 
in the theories. 

Note that we may use the bare potential Vq throughout 
these equations, since we consider the model Lieb-Liniger 
system in which the interaction potential truly is a delta 
function. Furthermore, in ID there is an absence of the 
ultra-violet divergences which occur in other dimensions 
when a contact potential is used. Since the above equa- 
tions are always solved self-consistently, we still generate 
the many-body T-matrix (where used) to all orders from 
the bare potential [T^ |. 

A. Condensate and quasi-condensate properties 

We start our simulations in the weakly-interacting 
limit with a small value of Vq. In this limit Lieb and 
Liniger 4] found that the Bogoliubov theory gave a good 
description of the homogeneous system, and so we expect 
our model should also be valid. We solve the above equa- 
tions self-consistently to obtain the wave function, chem- 
ical potential, many-body T-matrix, p(x), and k(x). We 
then ramp up the value of Vo, using the results from 
the previous iteration as our starting point. By this 
method we have managed to solve the equations from 
the weakly-interacting limit into the strongly-interacting 



limit where the system should become a Tonks-Girardeau 
gas. Again, we expect the chemical potential to change 

2/3 

initially as p oc V at low Vo , before reaching a plateau 
at p = N iot hu> in the strongly- interacting limit. 

We present the results from the different forms of the 
theory below. In all cases we have used Nq = 15, and 
we start our simulations with Vq = Vo/tiu)£t Ta ,p = 0.2 
which is sufficiently weak an interaction that a solution 
is obtained if p and k are initially zero. 

In Fig. the chemical potential has been plotted as 
a function of the potential strength Vo, for the various 
different models, including the GPE results from the pre- 
vious section. Figure |3Jd also shows p/N tot the chemical 
potential per atom, in order to compare more easily to 
the expected high-Vo limit p = N^Tkj, Note that A^ot 
is the total population of the system, consisting of both 
the condensate and non-condensate populations. In the 
low-Vo limit it can be seen that all of the models provide 
very similar results, as should be expected since the gas is 
weakly-interacting in this limit. As Vo is increased, how- 
ever, the different models show quite different behaviour. 

In terms of the chemical potential, neither the BdG 
nor the HFB-Popov models give the required behaviour, 
instead p increases rapidly as Vq is increased. In neither 
model does the graph show a plateau in the high-Vo limit, 
and the chemical potential greatly exceeds the Tonks- 
Girardeau limit p — N tot hu) in the HFB-Popov case. 

2/3 

Instead the chemical potential roughly follows the V 
behaviour for all Vq. The differences between these two 
models is seen in Fig. |5p, in which the results differ due 
to the different total populations. In the BdG case the 
non-condensate population increases very rapidly, such 
that in the very high-Vo limit p/N to t is actually decreas- 
ing. In the HFB-Popov model, the interactions between 
the condensate and non-condensate slow the growth of 
the non-condensate population, but the chemical poten- 
tial still rises unchecked. 

The disagreement of these models with the exact re- 
sults is not at all surprising. In both models the anoma- 
lous average is completely neglected, and so the many- 
body T-matrix is not introduced at all. We have argued 
here, and in earlier papers |?3.ll6|. about the importance 
of many-body effects in interactions in low-dimensions, 
and so we expect that neglecting these effects will lead 
to an inadequate theory. The BdG and HFB-Popov re- 
sults presented here are similar to those obtained in the 
homogeneous gas by Lieb and Liniger 4] , who found that 
the standard Bogoliubov method failed at high interac- 
tion strengths, when 7 > 2. In our results we also see 
quite good agreement with the GPE of the previous sec- 
tion for 7tr a p ^5 1 an d the differences between the theories 
become much greater when 7t rap is increased above this 
point j^. Note that 7t rap ~ 2 corresponds to Vq ~ 10 in 
the figure, which is roughly the point at which the curves 
separate. 

The anomalous average, and therefore the many-body 
T-matrix, can first be included by solving the GHFB1 
equations. Figure shows that there is little difference 
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FIG. 5: a) The chemical potential as a function of the inter- 
action potential Vb = Vo/hu)£tta.p as calculated by the various 
different models. The solid lines represent the results from 
the GHFB2 theory (upper line) and the GPE results of the 
previous section (lower line) . The other results correspond to 
the GHFB1 theory (dashed), HFB-Popov theory (dash-dot) 
and the BdG equations (dotted), b) The chemical potential 
per atom in the system as a function of Vo- The exact results 
predict that fi/Ntatfow = 1 in the high Vo limit. The bro- 
ken lines correspond to the same curves as in part (a), while 
the solid lines swap - the lower curve representing the GHFB2 
results in this figure. The parameter 7t rap is not plotted, how- 
ever for comparison purposes Vo = 200 roughly corresponds 
to 7 trap ft! 80. 



between these results and those of the HFB-Popov the- 
ory. The chemical potential still easily exceeds the limit 
of the Tonks-Girardeau gas, and so the GHFB1 theory 
also fails rapidly as Vb increases. 

However, the GHFB1 results do allow us to examine 
the form of the many-body T-matrix. In Fig.[S]the spa- 
tial dependence of the many-body T-matrix is plotted as 
Vb is increased. In fact, the value being plotted is the T- 
matrix divided by the condensate density |-0(a;)| 2 . From 
the previous section we expect that the T-matrix be- 
comes proportional to the condensate density as Vo — ► oo, 
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FIG. 6: TM_-e,{x)/N \%l)(x)\ 2 for the many-body T-matrix cal- 
culated from the GHFB1 theory for various values of Vb. Note 
that the ratio is only plotted for the range over which the con- 
densate density is significant. The values of Vb used were 0.2, 
2.4, 6.5, 15 and 44, which correspond to the narrowest to 
the widest curves respectively. These values correspond to 



7tr. 



= 0.02, 0.46, 1.6, 4.3 and 16.5. 



as in Eq. Q14JI . Figure shows that as the interaction 
strength is increased this becomes an increasingly good 
approximation. The very good agreement with the form 
expected for large Vq is further evidence that the cou- 
pling parameter ghh{x) derived earlier is accurate in this 
limit. This verifies our approximation for the many-body 
T-matrix in terms of the off-shell two-body T-matrix nu- 
merically. 

That the GHFB1 theory fails in the high-Vb limit is 
not particularly surprising. In the low-Vb limit almost 
all the atoms are in the condensate. As the interaction 
strength is increased, however, we reach the regime of 
the quasicondensate discussed by Petrov et al. [l(| • As 
Vb is increased the number of atoms in excited states 
(i.e. p) increases, and can become a sizable proportion 
of the total number of atoms in the system. It therefore 
becomes important to describe the interactions between 
these atoms accurately. 

The description of interactions between excited atoms 
is more complicated than those between condensate 
atoms. However, two limits may be investigated rela- 
tively straightforwardly. If the majority of the collisions 
occur at a high energy then the correct description would 
be given by the two-body T-matrix (since the many-body 
T-matrix tends towards the two-body T-matrix in this 
limit). On the other hand, if the majority of the colli- 
sions are at low energy, then it is more appropriate to 
describe them by the same many-body T-matrix as used 
to treat the condensate-condensate interactions. In 2D 
and 3D it might be expected that the large density of 
states at high energies should mean that the majority of 
collisions are in this regime. However, in ID the density 
of states is independent of energy in a trap. In the zero- 
temperature limit in which we are working it therefore 
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seems probable that the appropriate description of the 
interactions involving non-condensate atoms is the same 
many-body T-matrix used to describe the condensate in- 
teractions. 

The many-body T-matrix is used to describe non- 
condensate interactions in the GHFB2 version of the the- 
ory. It can be seen from Fig. [5] that the results using 
this theory are quite different from those of the GHFB1 
theory. In the low-Vb limit the results are similar, but 
as the interactions are increased the chemical potential 
becomes independent of Vb, as we expect in the Tonks- 
Girardeau limit. The figure clearly shows that it reaches 
an asymptotic value of ~ 0.7Nt o thu). This is lower than 
the expected value of N tot hcu, but nonetheless the results 
are in qualitative agreement with the behaviour expected 
for the Tonks gas. 

Again, we can compare the spatial shape of the many- 
body T-matrix with our earlier approximation, and the 
results look essentially the same as those calculated from 
the GHFB1 results (Fig. |HJ. The major change in 
the GHFB2 theory is therefore to provide a better de- 
scription of the condensate-non-condensate interactions 
rather than any implicit change in the form of the many- 
body T-matrix. 

We can also compare the predicted density distribution 
in the high-Vb limit with the solutions of the quintic GPE 
which was shown to give good agreement with the Tonks- 
Girardeau results in section lTTTl In Fig. [7] we compare the 
total densities [-/Vol "0(2)1 2 + p(x)) for the different models 
with the Thomas- Fermi prediction from Eq. (|16fl . It can 
easily be seen that the best agreement is given by the 
GHFB2 theory, while the BdG, HFB-Popov and GHFB1 
results all give entirely the wrong shape. It must be said 
however, that the GHFB2 results are still not in good 
agreement with the predicted density. 

Thus the GHFB2 theory interpreted in the normal 
manner fails to quantitatively describe the ID system 
in the impenetrable limit. This is not surprising since 
it assumes the existence of a condensate while we know 
that a condensate does not exist in this limit. Never- 
theless, the GHFB2 theory does give qualitatively the 
correct behaviour for the chemical potential and density 
in the high-Vb limit, which is a significant improvement 
on the other forms of mean-field theory used for 3D con- 
densates. 

The failure of the GHFB1 theory and the relative suc- 
cess of the GHFB2 theory indicate that the description 
of the non-condensate is the key to solving the ID prob- 
lem in the strongly- interacting limit. The GHFB1 model 
fails at high Vo since an increase in Vo leads to an in- 
crease in p. The p term obviously becomes far too domi- 
nant in GHFB1 and therefore needs to be suppressed by 
some means. In GHFB2 the many-body T-matrix is in- 
troduced to describe interactions involving p, and since 
Tmb/Vq decreases with increasing V this achieves the 
required suppression. As confirmation of this argument, 
the GHFB1 results were rerun with p artificially set to 
zero. With this adjustment, the GHFB1 model gave re- 
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FIG. 7: Total density N \ip{x)\ 2 + p{x) (solid) compared 
to the Quintic GPE prediction in the Thomas-Fermi limit 
(dashed) in the high Vo regime for the different theories, a) 
GHFB1 results at Vo = 50 (7t rap = 20) (this was as high 
a potential for which consistent results were obtained), b) 
GHFB2 results at V = 100 (7tra P = 30). c) HFB-Popov 
results at Vo = 20 (7tra P = 5) (again, results for potentials 
higher than this were not obtained), d) The BdG equation 
results at Vo = 100 (7t rap = 15). 



suits which looked qualitatively similar to the GHFB2 
results, with the chemical potential reaching a plateau in 
the high-Vb limit. 



B. Quasiparticle spectrum 

Having seen qualitative agreement between the 
GHFB2 theory and the expected behaviour of the chem- 
ical potential and the density, we turn now to investigate 
the quasiparticle spectrum. In the high-energy limit the 
quasiparticles are predicted to be single-particle like, and 
therefore for a ID harmonic potential the energy will rise 
linearly with the quasiparticle quantum number. How- 
ever, for the low-energy quasiparticle states (with ener- 
gies less than ~ p), we can use the hydrodynamic theory 
to obtain an expression for the quasiparticle energies e n 
in the Thomas-Fermi limit |26j| . Using the same method 
as in [2(j , we find that the low-energy quasiparticle spec- 
trum for a BEC in ID described by a cubic GPE is 



e n = huj 



n(n + 1) 



(28) 



where n = 1, 2, ... and n = 1 is the Kohn mode. However, 
we have shown earlier that in the impenetrable limit the 
system is described instead by a quintic GPE. For such 
a GPE the same hydrodynamic approach instead pre- 
dicts m 



nhui, 



(29) 
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which is obviously what one would expect for non- 
interacting fermions, and so we would expect this rela- 
tionship from a correct theory in the high-Vo limit. 

Figure GO shows the GHFB2 results for the energies of 
the first four quasiparticle modes as a function of the in- 
teraction strength. In the weakly-interacting limit the 
energies are all close to the prediction of Eq. Q28[). as ex- 
pected. However, as the interaction strength is increased, 
the quasiparticle energies do not go over to the limit of 
Eq. J2HJ). Instead they all fall to around 60% of that 
value and, even in the high-Vo limit, the relationship be- 
tween the different energies scales more like Eq. (|28[1 than 
Eq. H29f> . The conclusion which must be drawn is that 
the GHFB2 theory does not predict the quasiparticle en- 
ergies at all accurately in the high-Vb limit. This can be 
seen particularly strongly in the case of the Kohn mode 
(ei), which should be equal to Tilo at all interaction po- 
tentials since it simply corresponds to a centre-of-mass 
oscillation of the entire system in the harmonic poten- 
tial. 

The reason that GHFB2 does not predict the quasipar- 
ticle spectrum well is because it does not include any dy- 
namics of the non-condensate. In the strongly-interacting 
limit a large proportion of the atoms in the system are in 
the non-condensate, which is assumed to be completely 
static in all of the theories discussed here. Since the Kohn 
mode corresponds to the motion of all of the atoms to- 
gether, it is not surprising that a model which does not 
include the dynamics of a large proportion of those atoms 
predicts the wrong frequency. 

To see this in more detail, all of the theories discussed 
here can be derived by starting from the generalised time- 
dependent GPE 
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and expanding ty(x,t) in terms of the unperturbed 
ground-state ip(x) and time-dependent fluctuations 
Aip(x,t), and linearising (see, for example, j2^|). The 
trouble occurs when it is recognised that in ID the ap- 
propriate coupling parameter U con (x) is the many-body 
T-matrix, as argued earlier. In the high-Vo limit, the 
many-body T-matrix is proportional to N \ip(x)\ 2 in the 
static case as we have shown. Applying the local density 
approximation, the time-dependent form of the many- 
body T-matrix should then be No\fy (x , t)\ 2 (assuming 
that the time-scale for interparticle interactions is much 
quicker than that for condensate dynamics). In the lin- 
earisation procedure, one should therefore also expand 
the many-body T-matrix as 



T M b{x) (x \if)(x) + AV(M)f 



(31) 
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FIG. 8: The energies of the first four quasiparticle states 
as a function of the interaction potential Vo, as predicted by 
GHFB2 theory. The dashed lines indicate the value expected 
from hydrodynamic theory by the linearisation of the cubic 
GPE (for the Kohn mode [si] the expected value is lhco). The 
correct values in the high-Vo limit are hu, 2hui, 3hu>, and ATiuj 
(from left to right). 



earlier and hence to the result of Eq. I|29|l . We see there- 
fore that including the dynamics of the anomalous aver- 
age in the GHFB2 theory will lead to the correct predic- 
tion for the excitation spectrum in the high-Vo limit. For 
a general value of Vo, including the fluctuations of the 
many-body T-matrix requires a dynamic treatment of 
the anomalous average and hence of the non-condensate 
generally. 



C. Coherences 

It is of interest to calculate the coherence functions 
for the Lieb-Liniger gas from the GHFB2 results. We 
have argued that we are dealing with a gas in which a 
condensate may form in the low-Vo limit, but where it 
disappears when the interactions are increased. Since 
coherence is a major feature of a condensate we should 
be able to see a loss of coherence as Vo is increased, and 
again we can try to relate this to the expected results in 
the Tonks limit. 



1. First- order coherence 

The first order coherence function (r, is defined 
in terms of the quantum field operators by |29j 



Making this replacement for U con (x) in Eq. (|30|1 while 
setting p = leads directly to the quintic GPE discussed 



G (1) (r,r') = (*t( r )$( r ')) ) 



(32) 



11 



and the normalised first order coherence function is there- 
fore 
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(i) 



(r,r') = 



G«(r,r') 



y/GW(r, r)GW(ry) 



(33) 



Decomposing the field operator in terms of the conden- 
sate and its fluctuations leads to an expression for the 
coherence function at T — in terms of the conden- 
sate wave function and Bogoliubov quasiparticle ampli- 
tudes [13 

G«(r,r') = Ar ^(r)Vo(r') + ^«i(rK(r'). (34) 

Figure shows the normalised first order coherence 
function obtained from the GHFB2 results for various 
values of Vq. Note that the figure has been plotted in 
terms of the distance scale x = x/ Rtf, where i?TF is the 
Thomas-Fermi radius of the condensate, in order to take 
into account the greater spatial extent of a condensate 
with stronger interactions. Figure |H1 clearly shows that 
the system is coherent over its entire length in the weakly- 
interacting limit. However, as the interaction strength is 
increased the coherence can be seen to diminish substan- 
tially, and the width of the central maximum decreases 
as well. A diminishing range of coherence is consistent 
with the absence of a true condensate in the strongly- 
interacting regime. However, the rate at which the coher- 
ence calculated here decreases is surprisingly slow. Gan- 
gardt and Shlyapnikov [31| have presented results for the 
coherence function obtained from hydrodynamic theory 
applied to a ID gas, and they see a much greater fall in 
the range of the coherence length - dropping by almost 
an order of magnitude as 7t ra p decreases from 0.1 to 10. 
However, their results were obtained for a system of 10 
atoms, many more than we deal with here. Also, even for 
the largest value of Vq plotted, the condensate depletion 
is still only about 50%, in our case. 

2. Second-order coherence 

We now look briefly at the second-order coherence 
function defined by 



5 (2) (r,r') 



(tft(r)tft( r ')*( r ')*( r )) 
(*t( r )#( r ))(VEft(r')#(r'))' 



(35) 



We will consider this function only in the equal-position 
limit (i.e. when r = r'), in which case it may be written 
as 

g®(x,x) = (36) 
2N \^(x)\ 2 [p(x) + k(x)} + l/^)! 2 + \ K (x)\ 2 
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[N \^x)\2 + P (x)Y 



following a similar argument as for the first-order coher- 
ence function. The second-order equal-position coher- 
ence function is of interest because of the various limits in 




x/Rq 



FIG. 9: The normalised first order coherence function 
g^'(0, x) calculated using the GHFB2 theory for various val- 
ues of Vb- The solid curves represent (from top to bottom) 
the results for V = 0.2, 2.4, 6.5, 15, 100 and 200 (correspond- 
ing to 7tra P = 0.02, 0.46, 1.6, 4.1, 31 and 61 respectively). The 
coherence function is plotted as a function of x/ Rtf where 
Rtf is the Thomas- Fermi radius of the condensate. 



which exact results are known. For a pure BEC with per- 
fect second-order coherence we have g^ (x, x) = 1, whilst 
for a thermal cloud of bosons g^ 2 \x, x) — 2 (2^, 13^ . 
In the strong-coupling limit in ID however, the Bosc- 
Fermi mapping theorem predicts that g^ (x, x) = fl3| 
for impenetrable bosons. We would therefore expect to 
find a second-order coherence of close to 1 in the weakly- 
interacting limit, which decreases to as the interaction 
strength is increased, in contrast to the higher dimen- 
sional case where the coherence function would increase 
to 2 as the non-condensate became significant |39| . 

Figure ITUI shows g^ 2 \x,x) as calculated from the re- 
sults of the GHFB2 theory. The results can be seen to 
be close to unity, as expected, for the weakly-interacting 
case (and rising to 2 outside the region where the conden- 
sate exists). As the interactions are increased, <?( 2 ' (x, x) 
initially decreases closer to the g^ (x, x) = limit ex- 
pected. However, in the very strongly-interacting limit 
the coherence function rises again to become greater than 
1. Thus the GHFB2 theory, whilst in qualitative agree- 
ment with the behaviour of the chemical potential and 
the density, does not give the expected results for the 
coherence of the system in the impenetrable limit. 



V. REINTERPRETING THE RESULTS 

As seen in the previous section, although the GHFB2 
theory seems to go further than other standard theories in 
describing the ID Bose gas into the strongly-interacting 
limit, there are significant differences between the results 
in the impenetrable limit and the exact results of the 
Tonks-Girardeau gas. However, we will see in this sec- 
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FIG. 10: Plot of the normalised equal-position second-order 
coherence function g^ 2 \x,x) obtained from the results of 
GHFB2 theory for various values of Vb- The dashed line 
shows the weakly interacting limit with Vb = 0.2, the dot- 
dashed line is for Vb = 2.4. The dotted lines then correspond 
to values of Vb = 6.5, 15, 100 from bottom to top, and the fi- 
nal solid line is in the strongly interacting limit with Vb = 200 
(these values correspond to 7t rap = 0.02, 0.46, 1.6, 4.1, 31 and 
61 respectively) . Note that the :r-axis is scaled in units of the 
Thomas-Fermi radius of the system. 



tion that the agreement can be improved significantly by 
reinterpreting the results in a different manner. The dis- 
agreement with the expected values for the asymptotic 
value of the chemical potential and with the shape of the 
density are both related to the non-condensate popula- 
tion in the form of p{x). Looking again at Fig.[S]one can 
see that the chemical potential in the GHFB2 case rises 
to an asymptotic value of approximately \bTioj. Although 
this result was obtained for a condensate population of 
-/Vo = 15, it is in excellent agreement with the expected 
result for a system with 15 atoms in total. By itself this 
may appear to be a co-incidence; however the shape of 
the density resulting merely from the condensate "wave 
function" ip(x) is plotted against the Thomas-Fermi ex- 
act prediction in Fig. 1111 and the two shapes can be seen 
to agree quite well. It seems that the best agreement 
would be midway between that in Fig. ^2 an d that in 
Fig. 03 (where the total density has been plotted). 

This seems to imply that the GHFB2 theory predicts 
the correct results in the strongly-interacting limit if one 
partly ignores the contributions from the non-condensate 
(these contributions are sizeable - the non-condensate 
population for the results plotted in Fig. ^2 is about 
No/2). This would have to be combined with a re- 
interpretation of the condensate density as the quasi- 
condensate density and, ultimately, as the total density. 
Some procedure of this kind is clearly required in ID be- 
cause the various theories discussed in the previous sec- 
tion are based on the existence of a condensate whereas 
in fact the condensate disappears as Vb increases. 



FIG. 11: Density due to the condensate population only 
(solid) resulting from the GHFB2 theory at Vo = 100 (7 tra p = 
30). The dashed line again shows the Thomas- Fermi predic- 
tion for the exact density for 15 atoms in the Tonks-Girardeau 
limit. 



As another illustration of the good agreement with ex- 
act results which can be obtained with such a reinterpre- 
tation, in Fig. ^] the second-order equal-position coher- 
ence function is plotted with results obtained from the 
GHFB1 theory with p(x) artificially set to zero (as men- 
tioned earlier, this approach gave qualitative agreement 
with the expected behaviour of the chemical potential). 
The results in this case clearly show the behaviour pre- 
dicted in the previous section, with g^ 2 \x,x) decreasing 
from unity to zero as the interaction strength is increased. 
Such a theory (with p(x) = 0) would result from a reinter- 
pretation of No\ip(x)\ 2 as the total density of the system, 
rather than merely that of the condensate. The agree- 
ment with the expected behaviour for (x, x) can easily 
be seen from the following argument. In the high-Vb limit 
we have shown that Tmb becomes independent of Vq. 
From Eq. I|27[l. this must mean that n(x) — > — NoiPq(x) 
in the same limit. Equation then clearly shows that 
the second order coherence function vanishes in the high- 
Vb limit if p = 0. 

That such a reinterpretation of the condensate is nec- 
essary and will lead to the correct mean-field theory at 
large Vb can be seen by comparison with the GPE results 
in section Till CI If we set p — and reinterpret na(x) as 
"tot (x) then the GHFB theory reduces to the earlier GPE 
theory, the only difference being that g^o (x) is calculated 
from Eq. (|27[) instead of Eq. (|19|l . Assuming the validity 
of the local density approximation, however, these two 
definitions will lead to the same results. We note that 
a theory has recently been proposed by Anderson et al. 
[33L l34| which in effect leads to a similar reinterpretation 
of the quantities p(x), k(x) and uq(x) by considering the 
effects of phase fluctuations to all orders, and this ap- 
proach may provide a microscopic justification for this 
procedure. 



13 







1 1 1 
1 1 

1 ' 

1 

1 / 


" — - _ _ 




1 , 

1 

1 / 






+■ - - - 7 

1 


\ 




1 






1 


'--.w 







x/Rtf 

FIG. 12: Plot of the normalised equal-position second-order 
coherence function g^ 2 \x,x) obtained from the results of 
GHFB1 theory with p(x) artificially set to zero everywhere, 
for various values of Vq. The dashed line shows the weakly 
interacting limit with Vo = 0.2, the dot-dashed line is for 
Vo — 3. The dotted lines then correspond to values of 
Vo = 7.5 and 32, while the final solid line is in the strongly 
interacting limit with Vo = 48 (these values correspond to 
Ttrap = 0.02, 0.65, 2.1, 11 and 17 respectively). Note that the 
re-axis is scaled in units of the Thomas-Fermi radius of the 
system. 



VI. CONCLUSIONS 

In this paper we have given a form for an approximate 
many-body T-matrix to describe the scattering which 
occurs in a zero temperature ID Bose gas. This ap- 
proximate T-matrix enables the ID Bose system to be 
described by a Gross-Pitaevskii equation which is easily 
solved. We have shown that the solutions to this equation 
for a trapped system are in excellent agreement with the 
density of a ID Lieb-Liniger model gas in the two limits 
in which exact results are known, and that the density 
distribution of such a gas changes markedly as the inter- 
actions become stronger. 

We have also attempted to push various higher-order 
theories of BECs from the weakly-interacting regime into 
the strongly-interacting regime for the Lieb-Liniger gas. 
We have shown that the theories which neglect many- 
body effects on the interactions (BdG and HFB-Popov) 
fail to provide adequate results at even quite low interac- 
tion strengths. The GHFB1 theory includes the many- 
body T-matrix for condensate interactions, but it also 
fails to describe the strongly-interacting limit since it 
does not account for the condensate-non-condensate in- 
teractions properly and hence the terms involving p be- 
come too large. 

However, the GHFB2 theory, which includes many- 
body effects in the condensate-non-condensate interac- 
tions, gave results which were qualitatively correct even 
in the strongly- interacting limit, predicting the correct 



behaviour of the chemical potential and density profile. 
Quantitative agreement was not found, but the results 
show that GHFB2 describes the strongly-interacting 
limit to a much better degree than the other mean-field 
theories tested. 

GHFB2 theory still fails to adequately describe some 
aspects of the system however. The first-order coherence 
function shows qualitatively the behaviour expected, but 
the second-order coherence, while initially showing the 
correct behaviour in the low- to mid-Vb range, deviates 
from the expected limit at high-Vb. Furthermore, the 
quasiparticle excitation spectrum does not go over to the 
correct behaviour at high-Vb- 

The quasiparticle spectrum can be corrected if the the- 
ory is modified to allow for fluctuations in the many-body 
T-matrix, as well as in the condensate density. This 
requires a dynamical treatment of the non-condensate. 
At the same time, it seems that the key to pushing the 
GHFB2 theory into the high-Vb limit (obtaining the cor- 
rect coherences and quantitative agreement in general) 
lies in a reinterpretation of the non-condensate density 
as the total density, and we have presented some heuris- 
tic evidence that this will lead to the expected results. 
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APPENDIX A: THE MANY-BODY T-MATRIX 
IN 3D 

As an illustration of our approximation of the many- 
body T-matrix in terms of the off-shell two-body T- 
matrix, we show in this appendix that in the three- 
dimensional case it reproduces the known results for the 
many-body T-matrix. 

The 3D case differs from ID and 2D because the many- 
body T-matrix is only a small correction to the two-body 
T-matrix. As mentioned earlier, the 3D two-body T- 
matrix is a constant Uq = Airh 2 a/m to first order in the 
low-energy limit. The first approximation to the 3D cou- 
pling parameter is therefore g^D = Uq. We can then ask 
if the method used in this paper for ID (and used in ref- 
erence |(| for 2D) accurately predicts the effects of the 
medium on the scattering. 

In 0] we showed that our approximation for the many- 
body T-matrix in 3D is 

16 

ff 3D = (0|T MB (£ = 0) 0) « (0\T 2B (E = 5 a* 0), 

(Al) 



14 



and the off-shell two-body T-matrix for a gas of hard- 
spheres of radius a was shown in |15| to be 



0\T™(E)\0) = 
Uo 
Uo 




1 mE 2 



for E > 0, 



1 mE 2 



(A2) 



for E < 0, 

Combining these equations gives (to first order in the 



small parameter k^a = dm^,a 2 /% ) 



ff3D = U 



( 16/xma 2x ' 



0[{k^af 



(A3) 



We can rewrite this in terms of the density no, since 
in a homogeneous system /i = g^T>no to leading order. 
Making this substitution and solving for g^o gives the 
result 



93D = Uo 



-.(n a 



3U/2 



O(n a 3 ), 



(A4) 



where we have expanded in terms of the dilute gas pa- 
rameter noa 3 , which is small for a dilute gas BEC. This 
result agrees with other solutions [25, 36] for the many- 
body T-matrix in 3D homogeneous systems. At T = we 
can, therefore, account for many-body effects on particle 
scattering using an off-shell two-body T-matrix. 
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